Indications of Spin-Charge Separation in the 2D Extended t-J Model 
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The 2D extended t — J model is studied computationally in a broad region of parameter space, mo- 
tivated by recent photoemission experiments for the undoped cuprate Ca2CuC>2Cl2 (F. Ronning et 
al., Science 282, 2067 (1998)). The one-hole ground state is shown to develop robust antiferromag- 
netic (AF) correlations between spins separated by the mobile hole (i.e. across the hole). This effect 
tends to decouple charge from spin, and the quasiparticle weight becomes negligible, particularly at 
momenta (0, n) — (tv, 0). Studies with more holes show precursors of metallic stripe formation, with 
holes sharing their individual spin arrangements, and AF-correlations generated across the stripe. 

PACS numbers: 74.20.-z, 74.20.Mn, 75.25.Dw 
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Electronic strong correlations are widely believed to 
be crucial for the explanation of the anomalous prop- 
erties of high-temperature superconductors. Several as- 
pects of the cuprate phenomenology are indeed contained 
in the two-dimensional (2D) t — J model, including ten- 
dencies towards d-wave pairing upon hole-doping. How- 
ever, pioneering angle-resolved photoemission (ARPES) 
studies by Wells et al. [Q on the insulating compound 
Sr2Cu02Cl2 revealed important discrepancies between 
experimental data and t — J model predictions Q] near 
momenta k = (tt,0) — (0, tt), of relevance for the doped 
cuprates. Moreover, recent ARPES experiments by Ron- 
ning et al. H for Ca2CuC>2Cl2 reported indications of a 
hole dispersion with d x 2_ y 2 -characteristics even in such 
an undoped compound, a remarkable and unexpected re- 
sult. Clearly the one-hole case must be better understood 
before addressing a finite hole density. 

In addition, 2D t — J model studies have not been able 
to stabilize the metallic stripes proposed as an explana- 
tion of neutron scattering experiments, with individual 
stripes affecting only one CuO chain with hole density 
n=0.5 Instead, a pattern with stripes involving two 
adjacent chains has emerged using the DMRG method 
with suitable boundary conditions J5| . Since the study of 
stripes in the 2D t — J model seems a subtle problem, it is 
important to develop new scenarios to generate metallic 
stripes upon doping to guide the interpretation of ex- 
perimental results. With a similar motivation, Laughlin 
recently argued for the possible existence of a new, 
but still unknown, fixed-point in an extended parameter 
space that could influence on the behavior of holes in 
antiferromagnets. Observing a new fixed-point is poten- 
tially important for a proper analysis of the cuprates. 

To search for theories beyond those currently available, 
a systematic analysis of the "extended" t— J model should 
be carried out. In this model, next-nearest-neighbor 
(NNN) hopping terms at distances \f2a (t') and 2a (t") 
are added to the standard t — J model, which only con- 
tains a nearest-neighbor amplitude t (a is the lattice 
spacing). The importance of NNN-hoppings to repro- 



duce ARPES results was discussed by Nazarenko et al. 
0, Belinicher et al. §, Eder et al. |g, Kim et al. [jioj, 
and addressed by other groups O). Tohyama et al. also 
remarked the importance of NNN-amplitudcs [^2) . The 
optimal values of t'/t,t"/t (-0.35 and 0.25, respectively 
H ) are compatible with band structure calculations [jl3| . 
It is currently accepted that the extended t — J model pro- 
duces a quasiparticle dispersion in excellent agreement 
with ARPES data. However, intuition is still lacking on 
the effect of t', t" on the behavior of holes. 

The main motivation of this paper is to contribute 
to the clarification of the physics contained in the 2D 
extended t — J model at low hole-density. Pursuing 
such goal, several surprises have been found. The most 
interesting is the stabilization by NNN-hoppings of a 
dynamically generated complex structure around mo- 
bile holes containing robust AF-correlations "across-the- 
hole" (AH) (Fig. la). Similar correlations were noticed in 
the standard t — J model on ladders [|l4| and using small 
clusters with NNN-hoppings |l5|. However, the physi- 
cal origin of the AH-correlations, the full spin arrange- 
ment around the hole, and specially its consequences, 
remain to be identified. In particular, due to the AH- 
structure here it is argued that the vicinity of the hole 
carries a small spin, correlated with a tiny quasiparti- 
cle (q.p.) weight Z in the one-particle spectral function, 
particularly at k = (it, 0) — (0, 7r). This is suggestive 
of spin-charge separation in two-dimensions, at least at 
short distances, a phenomenon searched for since the 
early proposals of high-T c theories. The present results 
suggest that this effect may be relevant at the couplings 
currently believed to be realistic for the cuprates. 

Our analysis starts with the observation that the one- 
hole ground-state q.p. weight Z in the t — J model van- 
ishes only in the limit when J/t— >0 since the size of 
the spin-1/2 spin-polaron around the hole, regulated by 
string excitations, diverges as J — * ^j. However, even 
at J/t=0 a state with vanishing Z is not obtained since 
in this limit the one-hole ground-state becomes ferromag- 
netic (FM) due to the Nagaoka mechanism. However, ad- 
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ditional hole-hopping terms may stabilize a more exotic 
ground-state since the extra hole-mobility can scramble 
severely the spin background near the hole, avoiding the 
localizing tendencies of the string excitations. 
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FIG. 1. (a) Schematic representation of the AF-bonds 
across the hole (solid lines) of the extended t— J model. Shown 
are five sites, with the hole at the center; (b) Lines of con- 
stant Z corresponding to momentum (n,0), obtained using a 
20-sites cluster with one-hole, a grid of 16 points in parameter 
space, and smooth interpolations among them, providing suf- 
ficient accuracy for our qualitative discussion, t'/t" is fixed to 
-1.4. Points A, B, C and D are mentioned in the text (B corre- 
sponds to the currently accepted coupling for the cuprates). 
In the inset, Z for an 18-site cluster and k = (7r/3, 7r/3) is 
shown vs. |t'/t|, at J/t = 0.125; (c) Energy of a hole on 
a 20-sites cluster for the available momenta, relative to the 
one-hole ground-state energy, t'/t — —0.35, t"/t = +0.25 
and J/t's are indicated. The momenta are in units of 7r/5. 

To analyze such a possibility, Z has been calculated in 
the extended t — J model. Results are shown in Fig. lb, 
using Exact Diagonalization (ED) techniques on 18- and 
20-sites clusters ||]. The ratio t'/t" was fixed to —1.4 in 
most of our study, as suggested by fits of ARPES dis- 
persions for Sr2Cu02Cl2 |J. From Fig. lb it is clear that 
adding amplitudes t'/t < 0, t"/t > drastically reduces 
Z at k = (it, 0), making it virtually negligible in large re- 
gions of parameter space. The inset of Fig. lb shows that 
similar results are obtained close to k = (tt/2, 7r/2), al- 
though the reduction of Z with increasing NNN-hoppings 
is much stronger at (tt, 0). The shape of the constant - 
Z lines seems mainly regulated by the strength of t',t" 
relative to J, reasonable since t is renormalized to J by 
AF-fluctuations while t',t" are only partially renor- 
malized since they connect same sublattice sites. Note 
that the region of very small J/t of the t — J model, hid- 
den by the FM instability, is substantially expanded by 
the addition of NNN-hoppings. 

The remarkable isotropy around (tt/2, tt/2) found in 
ARPES experiments [Q,|| for Sr 2 Cu0 2 Cl2 is reproduced 
in the extended t — J model p[-[Il]j, and Fig.lc for a 20- 
sites lattice confirms that the four available momenta 
near (tt/2, tt/2) indeed have lower energy than (tt,Q), 
once the proper set of NNN-hoppings M is used. How- 



ever, Fig.lc also shows that similar results are found at 
smaller couplings J/t=0.05, keeping the NNN-hoppings 
the same Q. This suggests that a quasi- isotropic dis- 
persion around (7r/2, 7r/2) may exist in the broad param- 
eter region with small Z. In this respect, points B, C and 
D of Fig. lb share similar properties, different from the 
traditionally studied regime of point A. This will be a 
recurrent conclusion of the results analyzed below. 
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FIG. 2. Spin-correlations in the mobile hole reference frame 
obtained using a 20-sites cluster with one-hole, and k = (it, 0). 
(a) corresponds to J/t = 0.4 and t' = t" = 0, while (b) 
is for J/t = 0.2, t'/t = -0.35, t"/t = +0.25. The width 
of the lines is proportional to the strength of the AF-bonds. 
Only AF-correlations at distances a and 2a are shown; (c) 
Strength of the AF-correlation across-the-hole Sah for the 
ground-state of one- hole on a 20-sites cluster with k = (tt, 0). 
A positive (negative) result represents an AF (FM) correla- 
tion. Shown are data in the x-direction (open) and y-direction 
(full). The points A, B, C, D (Fig.lb) and E correspond 
to (J/t = 0.4, t'/t = 0.0,t"/t = 0.0), (0.4, -0.35, +0.25), 
(0.2, -0.35, +0.25), (0.05, -0.35, +0.25), and (0.05, -1,+1), 
respectively; (d) Same as (c) but on a 18-sites cluster with 
k = (7r/3, 7t/3). Here x- and y-directions are equivalent. 

To understand the drastic reduction of Z after adding 
NNN-hoppings, consider in Figs.2a-b the spin correla- 
tions in the hole reference frame. Shown are results for 
a 20-sites cluster exactly solved and momentum (tt, 0), 
of importance for the cuprates. Results are presented 
for points A and C of Fig.lb. Similar results have been 
obtained using the ORBA technique |l7| on clusters of 
size 4x8 with 1.5 x 10 6 states, and also with ED on 16- 
sites clusters, suggesting that finite-size effects are small. 
Fig. 2a shows that in the absence of NNN-hoppings the 
spin correlations present simple AF tendencies E3| . The 
hole carries a spin cloud and Z is finite Q. However, 
including extra hoppings the spin correlations are qual- 
itatively different (Fig. 2b). Note the presence of AF- 
correlations across the hole, in both directions. The hole 
motion induces a coupling between sites of the same sub- 
lattice. It is important to note that the across-the-hole 
AF-bond is not isolated but it is supplemented in both 
directions by other robust AF-bonds resembling dynam- 



2 



ically generated ID Heisenberg segments along each axis 
(Fig. 2b), individually weakly coupled to the rest of the 
spins. Figs. 2c, d show the strength of the AH-bond at sev- 
eral points in parameter space. Comparing with Fig. lb, 
it is apparent that the stronger those correlations are, 
the weaker Z is. 

The robust AF-correlations among pairs of spins near 
the hole suggest that the total spin-1/2 of the one-hole 
problem, defined on an even number of sites cluster, may 
not be located near the hole. This can be analyzed by 
calculating the local spin (S?) at site i in the mobile hole 
reference frame, with the overall constraint ^;(Sf) = 
1/2. At J/t=0.4, t'=t"=0, Fig.3a shows (Sf> around the 
mobile hole for k = (tt, 0). Here the spin distribution is 
nontrivial and some results are even negative along the 
direction of movement of the hole. However, the large 
spin next to the hole along the y-axis suggests that spin 
and charge are bounded. On the other hand, at J/t = 
0.1, t' = — t" = — t (Fig. 3b), here shown as an extreme 
but illustrative example, (Sf ) is considerably more spread 
than in Fig. 3a. A robust AH-spin arrangement near the 
hole appears correlated with local spin-charge separation 
tendencies, which is intuitively reasonable. 
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FIG. 3. (a) (Sf) around a mobile hole (open circle) on a 
4x4 cluster with one hole, total spin z-projection +1/2, 
J/t = 0.4, t' = t" = 0, and k = (tt,0). The area of the 
circles is proportional to (Sf). The gray circles denote a neg- 
ative (Sf); (b) Same as (a) but for J/t = 0.1, t' = -t" = -t; 
(c) Schematic representation of the 3- and 9-site clusters men- 
tioned in the text, the latter at J/t = 0.1, t'/t = —0.35 and 
t"/t = 0.25. The hole is at the center and the solid lines rep- 
resent the ground-state dominant AF-bonds, which are stable 
in a broad region of parameter space. 

To gain insight into the AH-state, consider first just 
3 sites in an open chain, containing one hole, one spin- 
up and one down (6 states) . Since J is not crucial in this 
analysis simply use J=0, and allow for the hole to move at 
distance a (2a) with amplitude t (t"). Two eigenstates 
are in competition for the ground-state. One has FM- 
correlations (Nagaoka state) and the lowest energy for 
t" < 0. However, if t" > a s pin-singlet state with 
energy E = |(-t" - yj{t") 2 + 8t 2 ) is stabilized (Fig.3c). 
The spin-singlet nature of this state leads in a natural 
way to AH-bonds. Then, the extra t' < and t" > 
hoppings have the important role of favoring the singlet 



state over the competing FM-state 19 20[|. 

Further semi-quantitative insight can be gained from 
the analysis of a 9-site cluster with open boundary condi- 
tions (Fig. 3c). Here there are 630 states in the zero total- 
S z subspace, and characteristics already similar to those 
found in the numerical analysis of Figs.2a-b were ob- 
served in its solution. The ground state of this small clus- 
ter is dominated by a hole at the center, in a large region 
of parameter space. The similarities of the ground-state 
correlations on 9-sites and on larger clusters suggest, once 
again, that the effects discussed here are mainly local in 
space and they can be observed in simple toy models. 

Qualitative understanding of the AH-state formation 
can be obtained from the lD-Hubbard model at U/t = oo 
plf . In this limit holes and spins fully decouple and the 
wave function in the spin sector corresponds to that of 
a lD-Heisenberg chain involving all spins, as if the holes 
were absent from the problem. In other words, a state 

such as | ... H 1 — OH h — K . . ) , with AF-bonds across 

the holes, has an important weight in the ground state. 
Explicit calculations performed as part of this effort have 
indeed shown that the ID t — J — t' — t" model presents 
AH-bonds quite similar to those in Fig. 2b. However, if 
a similar procedure is attempted in 2D, i.e. the effective 
removal of the hole from the system by linking the two 
vertical and horizontal neighboring spins with a Heisen- 
berg coupling J, frustration cannot be avoided. Then, a 
straightforward generalization to 2D of the ID results is 
expected to fail. However, the numerical results shown 
here lead us to suggest that a compromise can be found 
between the tendencies to produce a charge-spin decou- 
pling and the frustration that prevents it. This can be 
achieved by creating ID-like spin arrangements in both 
directions, as shown in Fig. 2b. Each of these dynami- 
cally generated "chains" individually resemble the results 
found in the ID U/t = oo Hubbard model Q. By this 
procedure, spin and charge can locally decouple in 2D, 
as suggested by the computational results. 

The spin arrangement around the 2D hole at k = (tt, 0) 
reported here is different from previous alternatives dis- 
cussed in the literature. Although averaging spin correla- 
tions around the hole would only indicate a reduction of 
antifcrromagnetism in favor of a "spin-liquid" state, the 
spin arrangement in this region has more structure than 
a mere bubble of weakly interacting spins would have. In 
addition, while naively it may appear that the AH-bond 
could be described as a spin-singlet formation in a short- 
range resonant-valence-bond (RVB) state, Fig. 2b shows 
that it is better to represent it as an AF-bond part of a 
ID Heisenberg segment, an unexpected result in a 2D sys- 
tem. Note also that moving holes using NNN-hoppings 
in a RVB-state would tend to leave the "length" of the 
short spin-singlets at just a, an effect opposite to what 
was reported here, namely NNN-hoppings were found to 
enhance the AH-bonds. Then, the AH-state represents a 
new paradigm for the visualization of the spin distortion 
around a hole in an AF background. 
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FIG. 4. (a) AF-correlations corresponding to the 4x4 
two-hole ground-state with PBC and J/t=0.4, t'/t = -0.35, 
t"/t = +0.25, for the case when the mobile holes (open cir- 
cles) are projected to be at distance 2a along the x-axis. 
The bond width is proportional to the strength of the 
AF-correlation. The horizontal AH-bond strength is close 
to a perfect singlet, but this is due to the boundary con- 
ditions on the 16-sites cluster. In the bulk, its strength is 
expected to be reduced roughly by a factor 2; (b) Same as 
(a) but with the holes at distance a; (c) Same as (a) but 
with the holes at a larger distance; (d) N(q) for a two-holes 
4x4 cluster. Squares, circles and triangles are results 
for (J/t = 0.4,t'/t = 0,t"/t = 0), (0.4,-0.35,0.25), and 
(0.1, -0.35, 0.25), respectively. 

In Figs.4a-c results are shown for two holes at point B 
of Fig. lb, the realistic couplings for the cuprates. Fig. 4a 
contains the spin arrangement for two holes on a 4 x 4 
cluster with PBC and at distance 2a, obtained from the 
zero-momentum ground-state using a suitable projector 
operator. With one hole at the origin, the most likely 
position for the second hole is precisely at distance 2a, 
unlike for t',t" = which is dominated by distance a. 
The row where the holes are located is clearly differ- 
ent from the rest, resembling a 4 x 1 cluster with two 
holes. This suggests that holes tend to move in a ID- 
path with density 0.5. These ID-paths resemble metallic 
stripes, conceptually different from the insulating stripes 
with one- hole per site found near the large J/t phase 
separated region of the standard t — J model when 1/r- 
Coulomb interactions are added |22]| . Fig. 4b is a nat- 
ural consequence of Fig. 4a, when holes move using the 
hopping t. An interesting property of Fig.4a-b is the ar- 
rangement of the spins not on the 4x1 stripe. They 
form robust AF-bonds across the path of the two holes, 
as observed experimentally p| and in ladder studies [|l4| . 
This effect is not present for holes with low-mobility near 
phase separation, but in the context discussed here it is 
natural since it evolves from the AH-bonds of individual 
holes. Fig. 4c corresponds to holes not along a ID-path, 
carrying the surrounding environments of isolated holes. 
Fig.4d contains the Fourier-transform of the charge cor- 
relations N(q), enhanced at (V, 0) due to configurations 
Figs. 4a, b. Calculations for four holes on 16-sites show 



patterns similar to those in Fig. 4. In particular, the con- 
figuration with holes forming a square of side 2a has a 
substantial weight in the ground-state and contains two 
of the ID-like paths of Fig. 4a, running in each direction. 

It is important to remark that for the standard t — J 
model with J/t between 0.1 and 0.5 and t' = t" = the 
correlations across-the-stripe are actually FM on the 16- 
sites cluster, although the spin correlation among the two 
spins of the two-holes 4x1 row is still AF and robust. The 
NNN-hoppings are needed to stabilize a structure with 
characteristics similar to those observed in experiments 
0. Actually, results quantitatively similar to Figs.4a-c 
were obtained also at point C of Fig. lb, and even in an 
extreme case such as J/t = 0.1, t' = — t, t" = t. 

Summarizing, here it has been reported that a mobile 
hole in the extended t — J model generates a complex spin 
arrangement in its vicinity containing AF-bonds across 
the hole. This is correlated with a small Z and indications 
of spin-charge separation at k = (jr, 0). Further work 
will clarify the relevance of the one-hole AH-bonds at 
finite hole density. However, studies of two and four holes 
on small clusters have already provided a glimpse of a 
possible new mechanism for metallic stripe formation, 
which does not rely on phase separation tendencies and 
it is operative at small J/t. 
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